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Abstract 

The persistent current of correlated electrons in a continuous one-dimensional ring with a single scatterer is calculated by 
solving the many-body Schrodinger equation for several tens of electrons interacting via the electron-electron (e-e) interaction 
of finite range. The problem is solved by the configuration-interaction (CI) and diffusion Monte Carlo (DMC) methods. 
The CI and DMC results are in good agreement. In both cases, the persistent current / as a function of the ring length L 
exhibits the asymptotic dependence / oc L~^~°' typical of the Luttinger liquid, where the power a depends only on the e-e 
interaction. The numerical values of a agree with the known formula of the renormalisation-group theory. 
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A one- dimensional (ID) interacting electron system 
cannot be viewed as a Fermi liquid of non-interacting 
quasi-particles. Away from the charge-density- wave in- 
stability, the interacting ID system is a Luttinger liq- 
uid in which the elementary excitations are of collec- 
tive Bosonic nature (see the review [1]). The Luttinger 
liquid properties affect the ID electron transport. An 
interesting example of such effect is the persistent elec- 
tron current in a mesoscopic ID ring. Magnetic flux </> 
piercing the mesoscopic ring gives rise to the persistent 
current which can be expressed (at T — OK) as [2] 

/ = -dEo{4>)/d<i,, (1) 

where Eq is the eigenenergy of the many-body ground- 
state. If the ring is clean and the single-particle dis- 
persion law is parabolic, the electron-electron (e-e) in- 
teraction does not affect the persistent current owing 
to the Galilean invariance of the problem. However, if 
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a single scatterer (an impurity or a weak link) is in- 
troduced into the ring, the non-interacting and inter- 
acting result are expected to differ fundamentally. For 
non-interacting spinless electrons in the ID ring con- 
taining a single scatterer with transmission probability 
l^fcF 1^ ^ li the resulting persistent current depends on 
the magnetic flux and ring length (L) as [3] 

I = {evF/2L) 14^1 sin(<j!.'), (2) 

where </>' = 2iT(j)/4'o, 4>o ~ h/e is the flux quantum, kp 
is the Fermi wave vector, and uf is the Fermi velocity. 
For a spinless Luttinger liquid the persistent current 
follows the power law / oc L~"~^ . More precisely, [3] 

J (X L~°~"^ sin((/>'), (3) 

where the power a depends only on the e-e interac- 
tion, not on the properties of the scatterer. The for- 
mula (3) can also be obtained [3] heuristically as fol- 
lows. Matveev et al. [4] analyzed, how the e-e inter- 
action renormalizes the bare transmission amplitude 
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ikp of a single scatterer in the middle of the ID wire. 

They derived the rcnormalized amplitude tk^, by using 
the renormalization-group (RG) approach suitable for 
a weakly-interacting electron gas (a <C 1). In the limit 
|i/cir|'^ *C 1, their result can be expressed in the form 

U^-hAd/Lf, (4) 

where L is the wire length, d is the spatial range of the 
pair e-e interaction V{x — x'), and the power a is given 
(for spinless electrons) by expression 

a = [1/(0) - V{2kF)]/2nnvF, (5) 

with V{q) being the Fourier transform of V{x — x'). If 
one replaces in equation (2) the bare ikp by the rcnor- 
malized amplitude (4), one recovers the Luttinger- 
liquid dependence (3). This approach is heuristic and 
single-particle, nevertheless, a is now expressed micro- 
scopically in the form (5), expected to hold for a ^ 1. 

In the Luttinger- liquid model, the physics of the low- 
energy excitations is mapped onto an effective field the- 
ory using Bosonization [3], where terms expected to 
be negligible at low energies are omitted. Within this 
model, the asymptotic dependence (3) was obtained by 
another approximation, using the analogy to the prob- 
lem of quantum coherence in dissipative environment 
[3]. It would thus be useful to avoid this approxima- 
tion as well as Bosonisation. This has been done in Ref. 
[5] , where the persistent current was calculated numer- 
ically by solving the ID lattice model with nearest- 
neighbor hopping and interaction. The formula (3) was 
confirmed for long chains and strong scatterers. 

In this work the persistent current of correlated elec- 
trons is studied in a continuous ID ring, not in the ID 
lattice. Using the configuration-interaction (CI) and 
diffusion Monte Carlo (DMC) methods, we solve the 
continuous many-body Sclirodinger equation for sev- 
eral tens of electrons interacting via the e-e interaction 

V{x — x') = Vb exp(— \x — x'\ /d). (6) 

Interaction (6) emulates screening and allows compar- 
ison with correlated models [3,4,5] which also assume 
the e-e interaction of finite range. The figures 1 and 2 
present our major CI and DMC results. These results 
show the persistent current exhibiting the Luttinger- 
liquid law I oc L~^~" , already verified by the lattice 
model [5]. However, it has sofar not been verified, that 
the power a obeys the RG formula (5) also in the mi- 
croscopic model avoiding the RG approach. The CI and 
DMC data in figures 1 and 2 verify this fact quite well. 
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Fig. 1. Persistent current LI{(f> = O.250o) in the GaAs 
ring with a single scatterer as a function of the ring length 
L. The upper horizontal axis shows the number of elec- 
trons N, the electron density is N/L = 5 x lO'^ m~^. 
The transmission of the scatterer is = 0.03. The 

range of the e-e interaction is d = 3nm, the magnitude 
Vo is varied. The CI and DMC data are shown by sym- 
bols. The full lines connecting the DMC data arc a guide 
for eye. The dashed lines show the Luttinger liquid asymp- 
totics LI (X L~", with a obtained from the RG for- 
mula (5): a = 0.0277, 0.0561, 0.085.5, 0.1710, and 0.2565 for 
Vb = 11,22,34,68, and 102 meV, respectively. 

We consider N electrons with free motion along a 
circular ID ring threaded by magnetic flux due to the 
constant magnetic field. The non-interacting single- 
electron states (fi (x) obey the Schrddinger equation 



d^ 



+ ^S{x) ) <fii{x) = eiifii{x), 



(7) 



2m dx'^ 

with boundary condition (^i (a; -|-L) = exp (i^^^ 'fiiix), 
where x is the electron coordinate along the ring, m is 
the electron effective mass, and ')5{x) is the potential 
barrier due to the scatterer. The wave functions ipi{x) 
and energy levels ei arc found numerically as in Ref. [6]. 
We want to solve the Schrddinger equation H'^ = E'^ 
for N interacting electrons with Hamiltonian 



h' d^ 



H=y , o . 2 

^ \ 2m dxt 



(8) 



+ 'ySixi)\ 

Consider first the non-interacting many-body problem 
= ^Xi where H is Hamiltonian (8) without the e-e 
interaction. Clearly, in this case £i = £ii H \-Un ^^'^ 
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Fig. 2. Persistent current L/(</> = 0.25</)o) in the GaAs ring 
with a single scatterer as a function of the ring size. The 
parameters of the calculation are given in the figure. The 
CI and DMC data are shown by symbols. The dashed lines 
show the decay LI oc L~°', with a given by the RG formula 
(5): a = 0.2565 for d = 3nm and Vq = 102meV as well as 
for d = 12nm and Vb = 12.8meV. The CI and DMC data 
confirm this a. They also confirm that a remains unchanged 



when the transmission \ti. 



of the scatterer is varied. 



where i is the quantum number labeling a specific set of 
A*' occupied levels e^^ , . . . , ei„ . The CI method [7] relies 
on expansion 'I' = coXo + cixi + C2X2 + • . • • Using this 
expansion and equation ~ {xi\E\'^) we get 

00 

J2(S3S^J+V^J)CJ=EC^, i = 0, 1 , . . . , OO , (10) 

3 = 



where Vij = 



KXil Y.'^i^k - ^t)\Xj)- We reduce the 
infinite number of single-energy levels to the finite 
one by introducing a proper upper energy cutoff. This 
reduces the infinite number of equations (10) to a cer- 
tain finite number M + 1. We get the finite system 
E^lo (^i-^y + cj = Ec„ where i = 0,1,..., M. 
This system determines the eigenvalues Ei and eigen- 
vectors (co, Ci, . . . , Cm) for I = 0,1, ... , M. To obtain 
the persistent current (1), we need only the ground- 
state energy iJi=o as a function of 0. We have obtained 
Eo by solving the system with program ARPACK. 

We have also calculated Eo by the DMC method [7] . 
Technical details of this calculation are the same as in 
Ref. [8], the persistent current is obtained from (1). 

Our major CI and DMC results are shown in figures 
1 and 2. These figures show in log scale the persistent 
current I//(0 — 0.25(j)o) as a function of the ring length 
for a GaAs ring with one strong scatterer. In log scale 
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Fig. 3. Persistent current versus magnetic flux and energy 
of the many-body ground state for the GaAs ring with pa- 
rameters N = A and L = 0.08 fira. The transmission of the 
scatterer is set to p = 0.03. The parameters of the e-e 
interaction are Vq = 102 meV and d = 3 nm. The CI and 
DMC data are shown by open symbols, the lines connect- 
ing points are a guide for eye. The data for the persistent 
current follow the dependence /(0) <x sin(27r(/)/</>o) [3]. 

the decay LI oc L~" is linear, with slope —a. The CI 
and DMC data in figure 1 confirm the expected [3,5] 
trend: The stronger the e-e interaction the shorter the 
system necessary to reach the LI oc L~" asymptotics. 
In figure 2, where only the strongest e-e interaction is 
considered, the decay LI oc L"" exists already in rings 
with four electrons. To see this Luttinger-liquid fea- 
ture in such few-electron system is a surprising result. 
Further, the CI and DMC data in figures 1 and 2 con- 
firm that a obeys the RG formula (5). Finally, figure 3 
shows, that the persistent current is a sine function of 
the magnetic fiux also in the CI and DMC models . 
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